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Abstract 

We present our analysis on the clustering properties of star-forming galaxies se- 
lected by narrow-band excesses in the Subaru Deep Field. Specifically we focus on 
Ha emitting galaxies at z = 0.24 and z = 0.40 in the same field, to investigate possible 
evolutionary signatures of clustering properties of star-forming galaxies. Based on the 
analysis on 228 Ha emitting galaxies with 39.8 < logI/(Ha) < 40.8 at z = 0.40, we find 
that their two-point correlation function is estimated as ^= (r/1.62to;5oMpc)~^'^^^°'°^. 
This is similar to that of Ha emitting galaxies in the same Ha luminosity range at 
z = 0.24, = (r/1.88lo;49Mpc)~^'^^^°°'^. These correlation lengths are smaller than 
those for the brighter galaxy sample studied by Meneux et al. (2006) in the same 
redshift range. The evolution of correlation length between z = 0.24 and z = 0.40 is 
interpreted by the gravitational growth of the dark matter halos. 

Key words: galaxies: distances and redshifts — galaxies: evolution — cosmology: 
large-scale structure of universe 



1. INTRODUCTION 

Formation and evolution of galaxies and its dependence on the environment, e.g., the for- 
mation of large-scale structure, are fundamental problems of the modern cosmology. Studying 
the evolution of the clustering of galaxies is useful to understand the evolution of galaxies and 

* Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan. 
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that of large-scale structures. Galaxies are considered to form within dark matter halos. Owing 
to the WMAP observations (Spergel et al. 2007), the cosmological parameters are well deter- 
mined and the evolution of the clustering of dark matter halos is now well understood both 
analytically (Mo & White 1996; Sheth & Tormen 1999) and from N-body simulations (Jenkins 
et al. 1998; Kauffmann et al. 1999; Springel et al. 2005). However, physical relationship 
between galaxy populations and the dark matter halos is not clear, since galaxy formation 
processes include gas cooling, star formation, and feedback. 

At the present universe, galaxies selected from optical observations are considered to 
trace the distribution of underlying dark matter halos (Verde et al. 2002). However, the 
strength of the clustering depends on properties of galaxy populations: e.g., luminosity (Zehavi 
et al. 2005), color or spectral type (Zehavi et al. 2002), and stellar mass (Li et al. 2006). The 
dependence of clustering on the properties of galaxies gives constraint on theories of galaxy for- 
mation. To constrain the model of galaxy formation, it is useful to investigate such dependences 
in detail. In this paper, we concentrate the dependence of clustering on the luminosity of galax- 
ies. Norberget al. (2001) fit the relative bias 6/6* as a function of L/L*, 6/6* = 0.85 + 0.15L/L*, 
where 6* is the bias for L* galaxies and the bias is defined as the square root of the ratio of the 
galaxy and dark matter correlation functions, 6= (^g/^Divi)^^^- Using the SDSS data, Tegmark 
et al. (2004) gave another fitting function, 6/6* = 0.85 + 0.15L/L* -0.04(M-M*). The relative 
bias at the low-luminosity end is smaller than the Tegmark's relation in Zehavi et al. (2005). 
Therefore the relative bias of low-luminosity galaxies is still controversial. 

Recently, the dependence of clustering on the luminosity is studied up to z ~ 1.5 from 
the VIMOS-VLT Deep Survey (Marinoni et al. 2005; Meneux et al. 2006; Polio et al. 2006), 
the DEEP2 Galaxy Redshift Survey (Coil et al. 2006), and the Canada- France Legacy Survey 
(CFHTLS, McCracken et al. 2008). Brighter galaxies are more strongly clustered than low 
luminosity galaxies and the dependence of relative bias on the luminosity at z ~ 0.8 is in 
agreement with that in the local Universe (Marinoni et al. 2005). The observed decreasing of 
the clustering strength is consistent with simple gravitational growth picture (Meneux et al. 
2006). However, the evolution of the clustering of low- luminosity galaxies is not clear till now, 
since those objects are too faint for spectroscopic surveys. 

In this paper, we present the clustering properties of Ha emitters ai z = 0.24 and z = 0.40, 
which are low-luminosity active star-forming galaxies in the Subaru Deep Field. Their mean 
absolute magnitude is (Mb) ~ —17.0, which is 1 mag fainter than the previous observations 
(Polio et al. 2006; Meneux et al. 2006). We therefore study the evolution of galaxies with low 
luminosities. Using the narrowband imaging survey, we select star-forming galaxies within the 
restricted redshift to fainter magnitudes. Throughout this paper, magnitudes are given in the 
AB system. We adopt a fiat universe with ^matter = 0.3, Qa = 0.7, and Hq = 70 km s~^ Mpc~^. 
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2. PHOTOMETRIC CATALOG 



In this paper, we use the official photometric catalogs of the Subaru Deep Field (SDF) 
project that is a very deep optical imaging survey using the Suprime-Cam (Miyazaki et al. 
2002) on the 8.2 m Subaru Telescope (Kaifu et al. 2000; lye et al. 2004) at Mauna Kea 
Observatories. The SDF is located near the North Galactic Pole, centered at a(J2000) 
= 13''24™38^9 and (5(J2000) = +27°29'25".9. Details of the SDF project are given in 
Kashikawa et al.(2004). The SDF official photometric catalogs are obtained from the SDF 
web site (http://step.mtk.nao.ac.jp/sdf/project/). These official catalogs contain 5 broadband 
(B,V,Rc,i' , and z') and 2 narrowband {NB816 and A^i?921) photometric data. In this work, we 
use the i'-selected catalog with 2" diameter aperture photometry. The PSF size in this catalog is 
0.98 arcsec (Kashikawa et al. 2004). Since the Galactic extinction is not corrected in the magni- 
tudes in the official catalog, we applied the Galactic extinction correction of E{B — V)= 0.017 
(Schlegel et al. 1998). A photometric correction for each band is Ab = 0.067, Ay = 0.052, 
Ar^ = 0.043, Ai, = 0.033, A,> = 0.025, Anbsw = 0.030, and Anb92i = 0.024. 

The narrowband filters used in SDF are NB81Q = 8150 A, AA(FWHM) = 120 A] 
and NB921 [Ac = 9196 A, AA(FWHM) = 132 A]. These narrowband ffiter data are used for 
searching Ha emitters aX z = 0.24 and z = 0.40; note that the Ha luminosity function and the 
angular correlation function of Ha emitters ai z = 0.24 in the SDF has been already studied by 
our group (Morioka et al. 2008). The limiting magnitudes for a 3a detection on a 2" diameter 
aperture are Biim,3<7 = 28.45, \4im,3<7 = 27.74, i?ciim,3<7 = 27.80, il^^^^^ = 27 A3, z^^^g^ = 26.62, 
NB816iira,3a = 26.63, and AS921iim,3a = 26.54 . 

3. RESULTS 

3.1. Selection of Ha emitters 

First, we describe the selection method of Ha emitters a.t z = 0.40, which are basically 
similar to Fujita et al. (2003), Ly et al. (2007) and Shioya et al. (2008). For emission- 
line galaxies a.t z = 0.40, redshifted Ha emissions enter the NB921 band. We therefore select 
emission-line galaxies as NB921-excess objects. In order to select NB921-excess objects, we use 
z' band as off-band continuum. Taking account of the photometric error, we select NB921-excess 
objects using the following criteria: 

z' - NB921 > max[0.1, 3a{z' - NB921)], (1) 

where 

3a{z' - NB921) = -2.51og(l - v/(/3<x^b..J' + (/a^JV/^VBssi). (2) 

We show these criteria in Figure 1. In order to avoid the influence of saturation of brighter 
objects, we adopt another criterion of NB921 > 20. For galaxies with z' > z[^^^^^, we use the 
lower-limit value {z' — Ai?921)iow.iimit = -Zum so- " Ai?921 for our sample selection. We then select 
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the NB921-excess objects to the iVi?921faintest = 25.63 which is determined by A^i?921faintest + 
3(t{z' — NB921)NB92i=NB92if^i,^t^^^ = z'i:^rn 3cr- Ahhough the median of the z' — NB921 is shghtly 
different from (~ —0.05, see Ly et al. 2007, for this offset), we do not apply any correction 
for the NB921 and z' magnitudes. We then find 2039 sources that satisfy the above criteria. A 
narrowband survey of emission-hne galaxies potentially detects galaxies with different emission 
lines at different redshifts. The emission lines that can be detected in NB921 passband are Ha, 
H/3, [O III] AA4959,5007, [O ll] A3727, Lya, and so on. In order to distinguish Ha emitters 
at z = 0.40 from such emission-line objects, we investigate their broad-band colors comparing 
observed 2039 emitters with model spectral energy distributions (Coleman, Wu, & Weedman 
1980). In Figure 2, we show the B — vs. Rc — z' color-color diagram of 2039 sources and the 
loci of model galaxies. We select Ha emitters by using the following criteria, 

> l-08(i?c-^') + 0.55, (3) 

and 

5-i?c > 2.80(i?c-^')-0.84. (4) 

To investigate how our selection criteria suffer from contamination of galaxies at different 
redshifts, we plot the colors of galaxies with a spectroscopic redshift [specifically, 28 Ha emitters, 
60 [Olll] emitters, 24 H/9 emitters and 2 [On] emitters presented in Cowie et al. (2004), and 2 
Ha emitter, 22 [Olll] emitters and 4 [On] emitters presented in Ly et al. (2007)] in Figure 2. 
The contamination of [Olll], H/3 and [Oil] emitters into Ha is 1/112 (< 1%). The contamination 
of Ha emitters into [Olll], H/3 and [On] emitters is 3/30 (10%). These results are considered 
to justify our selection criteria. Finally, we select 356 Ha emitter candidates. 

We have already studied on Ha emitters aX z = 0.24 in the SDF (Morioka et al. 2008). 
Using the SDF official photometric catalog, we have selected Ha emitters a.t z = 0.24. Our 
selection criteria are as follows: (1) 20 < NB81Q < 26.1, (2) iz - NB81Q > max[0.1, 3(t(zz - 
NB816)], (3) {B-V)> l.Q{V - R^) - O.l L {B -V) > 3.1{V - R^) - 0.9, and (4) {B-V)> 
0.8{Rc - z') + 0.2 k{B-V)> 2.5{R^ - z') - 1.2, where iz continuum is defined as iz = 0.57 fi^ + 
OASfz'- In total, 258 Ha emitters are selected. 

3.2. Ha Luminosity of Ha emitters 

Since the clustering properties depend on the luminosity of galaxies, we have to know the 
luminosity range of our sample. In order to obtain the Ha luminosity for each source, we correct 
for the presence of [N ll] lines. Further, we also apply a mean internal extinction correction to 
each object. For these two corrections, we adopt the flux ratio of /(Ha) //([N ll]AA6548,6584) = 
2.3 (obtained by Kennicutt 1992; Gallego et al. 1997; used by Tresse & Maddox 1998; Yan et 
al. 1999; Iwamuro et al. 2000) and A^a = 1 (Gallego et al. 1995). We also apply a statistical 
correction (28% for NB921 and 21% for NB816; the average value of flux decrease due to 
the filter transmission) to the measured flux of each objects because the filter transmission 
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function is not square in shape (Fujita et al. 2003). The Ha flux is given from the observed 
flux, /obg(Ha+[N II]), by: 

/cor(Ha) = f^^,{}la + [N ll]) x ^(^f ^[j^ „j) x lO'-'^-- xC, (5) 

where C = 1.28 for NB921 and 1.21 for NB816. Finally, the Ha luminosity is given by L(Ha) = 
47r(iL/cor(Ha) where c^l is the luminosity distance at the redshift corresponding to the center of 
the filter passband: = 2.17 Gpc for z = 0.40 and 1.22 Gpc for z = 0.24. 

Figure 3 shows the Ha luminosity functions of our samples. The Ha luminosity function 
at z = 0.24 is considered to be complete between logL(Ha) = 39.4 and logL(Ha) = 40.8. On 
the other hand, that at z = 0.40 is incomplete for logL(Ha) < 39.8. To compare the clustering 
properties of Ha emitters without the luminosity effect, we made a subsample with the same 
luminosity range [39.8 < logL(Ha) < 40.8]. The numbers of Ha emitters whose luminosity 
range from logL(Ha) = 39.8 to logL(Ha) = 40.8 are 139 at ^ = 0.24 and 228 at z = 0.40. We 
also made another subsample with logL(Ha) > 40.8 for Ha emitters at z = 0.40 which contains 
126 Ha emitters. 

3.3. Spatial Distribution and Angular Two-Point Correlation Function 

Figure 4 shows the spatial distribution of our Ha emitters in the SDF. Left panel shows 
that at z = 0.24 and right panel shows that at z = 0.40. Objects in the range of 39.8 < logL(Ha) < 
40.8 are shown with larger filled circles. To discuss the clustering properties quantitatively, we 
derive the angular two-point correlation function (ACF), w{6), using the estimator defined by 
Landy & Szalay (1993), 

DDi9)-2DRie)+RR{e) 

^(^^ = mo) ' ^ ^ 

where DD{6), DR{6), and RR{6) are normalized numbers of galaxy-galaxy, galaxy-random, 
random-random pairs, respectively. The random sample consists of 100,000 sources with the 
same geometrical constraints as the galaxy sample. The formal error in w{6) is described by 



a^ = ,J[l + w{e)]/DD (7) 

(Hewett 1982). Because of the finite size of the survey, this estimate will be negative offset 
from the true w{6), which is called the integral constant. We calculate the integral constant, 
C = 0.02 using the following definition (Roche et al. 2002; see also Kovac et al. 2007): 

^- ERR • 

This value is small and we neglect it in further calculations. Figure 5 shows the ACF of 
whole sample of 356 Ha emitters in the SDF. The ACF is fit well by power law, w{d) = 
0.0047inoo9^"°'^^^°'°^- In this section, we use data points between 0.001° and 0.1° for the 
power- law fit. 
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It is useful to evaluate the correlation length ro of the two-point correlation function 
^(r) = (r/ro)^'^. The correlation length is derived from the ACF through Limber's equation 
(e.g., Peebles 1980). Assuming that the redshift distribution of Ha emitters is a top hat shape of 
z = 0.401 ±0.010, we obtain the correlation length of ro = 1.51 Mpc. The two-point correlation 
function for all Ha emitters is written as ^{r) = (r/1.51Mpc)~^'^^. We also evaluate the ACF and 
correlation function for both subsamples of Ha emitters ai z = 0.40. The ACF of luminous Ha 
emitters [logL(Ha) > 40.8] is fit by power law, w{e) =0m80tomll0''°-^^^°-^^ and the correlation 
function is ^(r) = {l/lMtlHUpc)-^-^^^^-^^ . The ACF of low-luminosity Ha emitters [39.8 < 
logL(Ha) < 40.8] is fit by power law, w{e) = 0.0078l[J;[]|]ie-°-s^^° °^ and correlation function is 
^(r) = (l/1.62lJ]|^Mpc)-i-S4±o-°^ For Ha emitters at z = 0.24 (39.8 < logL[Ha] < 40.8), the 
ACF is w{e) = OM3t°ofie-^-^^^°-^'^ and correlation function is ^{r) = (l/1.88in^Mpc)-i-S9±°-°^. 
For Ha emitters at z = 0.24, we assume that the redshift distribution of them is a top hat shape 
of ^ = 0.242 ± 0.009. We show these ACFs in Figure 6. 

4. DISCUSSION 

Figure 7 shows the correlation length function of redshift. For comparison, we 

also show those at 2 < 0.15 in the 2dFGRS (Norberg et al. 2001) and at 0.2 < 2; < 1.5 in 
the VIMOS-VLT Deep Survey (VVDS) (Meneux et al. 2006) together with our results of 
Ha emitters at z = 0.40 and 0.24. We also show the correlation length of Ha emitters at 
z = 0.24 in the COSMOS field (Shioya et al. 2008). The two-point correlation function of Ha 
emitters at ^ = 0.24 in the COSMOS field, ^ = (r/1.87i|]|jMpc)-i-*^s±°-°3, is very similar to that 
for our sample at 2 = 0.24, although we note that their average 5-band absolute magnitude, 
(Mb) = —18.5, is 1.5 mag brighter than our sample. 

First, the correlation lengths of our Ha emitters [39.8 < logL(Ha) < 40.8] at z = 0.24 and 
0.40 are smaller than those of early- and late-type galaxies at 0.2 < z < 0.5 derived by Meneux 
et al. (2006) from VVDS. Since the correlation length depends on the luminosity of galaxies, we 
need to make a fair comparison using absolute magnitudes in the same rest-frame wavelength. 
The mean S-band absolute magnitudes, (Mb), of our sample with 39.8 < logL(Ha) < 40.8 
are {Mb) = —16.9 and (Mb) = —17.0 for z = 0.40 and z = 0.24, respectively. On the other 
hand, the mean 5-band absolute magnitudes of Meneux's sample are (Mb) = —19.1 and —18.1 
for early-type and late-type galaxies, respectively. The smaller correlation length of our Ha 
emitters is interpreted by their small luminosity. 

We note that the clustering amplitude of our luminous Ha emitter sample at z = 0.40 
is weaker than those of the samples of Meneux et al. (2006). It may imply that the small 
clustering amplitude of our Ha emitter sample depends not only on the low luminosity but 
also the spectral type, e.g., the clustering amplitude of red galaxies is larger than that of blue 
galaxies. As we shown in Figure 2, the color of most of our sample is consistent with Scd - Irr 
type in Coleman et al. (1980), which is the same as type 3 & 4 in Meneux et al. (2006). The 
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clustering amplitude may depends on the equivalent widths of emission lines rather than the 
global SED. Our sample is biased toward galaxies with large emission-line equivalent width, 
which are considered to be in a very active phase of the star formation (e.g., Leitherer et al. 
1999). We note that there is no significant difference in the correlation length between our two 
samples although it would be expected that the brighter sample could have a larger correlation 
length than that for the sample with 39.8 < logL(Ha) < 40.8. One possibility is that our 
brighter sample galaxies are due to brightening by the current star formation activity. 

Next, we compare the correlation length of the Ha emitters a.t z = 0.40 with that of the 
Ha emitters ai z = 0.24. Although these are consistent within the error, the ratio of the tq at 
z = 0.40 to that at z = 0.24 (0.86) is also consistent with the expected value from the gravitational 
growth of dark matter halos. This fact implies that the Ha emitters at z = 0.40 exist in the 
same environment where the Ha emitters at z = 0.24 exist. If we assume the gravitational 
growth of dark matter halos to 2 ~ 0, the correlation length of this kind of galaxies evolve to 
tq ~ 2.3 Mpc. The ratio of the correlation length to that of the L* galaxies (M* = —20.5) 
is 0.36. This ratio is smaller than the prediction of Norberg's law. It is suggested that the 
clustering amplitude of low-luminosity galaxies become smaller for lower luminosity galaxies. 

We would like to thank the Subaru Telescope staff for their invaluable help. We would 
like to thank Nobunari Kashikawa who is the chief of SDF, and the other SDF members. This 
work has been financially supported in part by grants of JSPS (Nos. 15340059 and 17253001). 
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Table 1. Best-fit parameters of Angular Correlation Functions. 



Sample 




P 




z = 0.40 (all) 

z = 0.40 (log L(Ha) > 40.8) 

z = 0.40 (39.8 < logL(Ha) < 40.8) 

z = 0.24 (39.8 < logL(Ha) < 40.8) 


0.0047llJ:OS]ig 
0.0078l°:°gg 

0.013^°:°°^ 


0.94 ±0.04 
0.86±0.15 
0.84 ±0.08 
0.89 ±0.07 


1 CI +0.28 
^•O-'^-0.28 

l-86tj;g8 

1 ««+0.60 
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20 25 
NB921 

Fig. 1. Diagram between z' — NB921 vs. NB921. The gray horizontal sohd Une corresponds to 
z' - NB921 = 0.1. The gray vertical dotted hnes correspond to NB921 = 20 and 25.63 (see text). The gray 
dashed line shows the distribution of 3cr error. The gray dot-dashed line shows the 3a limiting magnitude 
of z'. 
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12 3 

Rc-z' 




Rc-z' 

Fig. 2. (Top) Color loci of model galaxies (CWW) from z = to z = 3 are shown with dotted lines: red, 
orange, green, and blue lines show the loci of E, Sbc, Scd, Irr galaxies, respectively. Colors oi z = 0.40, 0.85, 
0.89, and 1.47 (for Ha, [O ill], H/3, and [O ii] emitters, respectively) are shown with red, orange, green, 
and blue lines, respectively. Galaxies in the GOODS-N (Cowie et al. 2004) with redshifts corresponding 
to Ha emitters, [Oiii] emitters, H/3 emitters, and [On] emitters are shown as red, orange, green, and blue 
open squares, respectively. Galaxies in the SDF (Ly et al. 2007) with redshifts corresponding to Ha 
emitters, [Oiii] emitters, and [On] emitters are shown as red, orange, and blue filled squares, respectively. 
(Bottom) Diagram between B ~ Rc vs. Rc — z' for the 2039 emitter sources. Our 356 Ha emitters are 
isolated above the two gray solid lines. 
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40 42 

log L(Ha ) (erg s~^) 

Fig. 3. Ha luminosity functions at z = 0.24 (open squares) and at 2 = 0.40 (filled circles). Since the 
Ha luminosity function at z = 0.24 is considered to be complete at 39.8 < logi(Ha) < 40.8, we make 
subsamples of Ha emitters both at 2; = 0.24 and z = 0.40 within the range of 39.8 < logL(Ha) < 40.8. 
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Fig. 4. (Left) Spatial distribution of Ha emitters at z = 0.24. Large filled circles show objects in 
39.8 < log(Ha) < 40.8. (Right) Spatial distribution of Ha emitters at z = 0.40. Large filled circles show 
objects in 39.8 < log(Ha) < 40.8 and small ones show objects in log(Ha) > 40.8. Shadowed regions show 
the areas masked out for the detection. The field size of the SDF is 37' x 27'. 
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e (deg) 



Fig. 5. Angular two-point correlation function of all Ha emitter candidates (crosses), the luminous 
subsample (open boxes), and the faint subsample (open circles) at z = 0.40. Solid line, dashed line, 
and dotted line show the best-fit single power law for all Ha emitters, the luminous subsample, and the 
low-luminosity subsample, respectively. 
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Fig. 6. Angular two-point correlation function of Ha emitter candidates of 39.8 < logi(Ha) < 40.8 (filled 
circles). Solid line shows the best-fit single power law: w{6) — 0.00786'~°-®'*. For comparison, we also plot 
the angular two-point correlation function of Ha emitters at z = 0.24 of the same luminosity range (open 
circles). Dotted line shows the best-fit single power law: w{0) = O.O130~°'®^. 
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Fig. 7. Correlation length as a function of redshift. Our results are shown as a blue open square (the 
luminous subsample at z = 0.40), a blue open circles (the low-luminosity subsample at z = 0.40), and a 
blue filled circle (Ha emitters at z = 0.24). Blue triangle shows the correlation length for Ha emitters at 
z = 0.24 in the COSMOS field (Shioya et al. 2008). Red circles, green diamonds, and cyan crosses show 
early- type galaxies, late type galaxies, and starburst galaxies in Meneux et al. (2006). Filled black circles 
show the aX z < 0.15 for various Mb (Norberg et al. 2002). The dotted and dashed curves are tq values 
as predicted by the "e-modcl" (Groth & Peebles 1977) with e = 0.8 (corresponds to linear theory) and 
7 = 1.9. The dotted curves show the predicted gravitational growth of the dark matter for galaxies with 
Mb ^ —20.9 and Mb = —19.0 of Norberg ct al. (2001). The dashed curve show the predicted gravitational 
growth of the dark matter for our Ha emitters at z = 0.24. The number near each mark shows its absolute 
B magnitude, {Mb)- 
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